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Abstract 

We study the ground-state properties of mixtures of strongly interacting bosonic atoms in an 
optical lattice. Applying a mean-field approximation to the Hubbard model for Bose-Bose mixtures, 
we calculate the densities and superfluid order parameters for both species. Due to the repulsive 
interaction between the two species, the system exhibits phase separation. First, in the extreme 
limit of the zero-hopping case, we derive analytical expressions for the phase boundaries. In 
particular, we derive the conditions for phase separation in the Mott insulator phase. We find that 
the conditions for the phase separation depend on the on-site interactions as well as the occupation 
numbers. In particular, we show that the coexisting state appears by varying the on-site inter- 
species interaction. We also show the phase diagram of the finite hopping case. Second, we calculate 
the spatial density profile of 87 Rb- 41 K mixtures in the combined potential of a parabolic trap and 
an optical lattice using the local density approximation. We fixed the number of Rb and varied 
the number of K, and used the parameters estimated by experiments. We show that the phase 
separated 87 Rb- 41 K mixtures distribute like in a parabolic trap case. Furthermore, we find that 
phase separated mixtures distribute a nesting structure. 
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I. INTRODUCTION 



Systems of ultracold atoms in an optical lattice have attracted attention for studies of 
strongly correlated quantum matter. One of the most interesting phenomena exhibited by 
ultracold Bose atoms in an optical lattice is the superffuid (SF) to Mott insulator (MI) 



phase transition, which has been experimentally observed [l| . The properties of this system 
are well-captured by the Bose-Hubbard model In a shallow optical lattice, bosonic 

atoms are in an SF phase, wherein atoms have long-range phase coherence and number 
fluctuations. On the other hand, atoms are in an MI phase in a deep optical lattice, wherein 
the filling factor is fixed to an integer and phase coherence disappears. 

By the addition of a second atomic species, the system exhibits rich quantum phases. 



Some ex per imental groups are working on Fermi-Fermi 
Bose 
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llj mixtures in optical lattices. Recently, t 



6], Bose- Fermi 7H9], and Bose- 
re Florence group experimentally showed 



10[. The observed interference pattern 



87 Rb- 41 K mixtures trapped in an optical lattice 
revealed that the visibility of 87 Rb is lower than the pure 87 Rb case. 

The properties of Bose-Bose mixtures in an optical lattice have been studied theoretically 
using various approaches, such as by applying the Hubbard model using the perturbation 



theory 12[, the Gutzwiller approach 



13 



and quantum Monte Carlo simulations 



, the time evolving block decimation method 
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151 ] . In these previous studies, a supersolid phase 



and a counterflow superfluid phase have been predicted to exist, assuming that the two 
components have the same mass. On the other hand, it is important to consider the difference 
in components to investigate experimentally obtained two-species mixtures in an optical 
lattice. 

In this study, we have investigated two-species Bose mixtures in an optical lattice. In 
particular, we focus on the conditions for phase-separation. We have considered realistic 
parameter sets for the 87 Rb- 41 K mixture. Specifically, we consider mixtures of different 
bosonic atoms with different masses, on-site interactions, and hopping amplitudes. Our 
main purpose is to clarify superfluid order parameters in a combined potential of a parabolic 
potential and an optical lattice. We apply the mean-field approximation and include the 
effects of trap potential by the local density approximation. We calculate spatial profiles of 
the densities and the superfluid order parameters. 

This paper is organized as follows: In section [III we explain the mean- field approximation 
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and the numerical method used. In section Hilt we present the conditions for the spatial 
phase separation of Bose-Bose mixtures in an optical lattice. In order to gain qualitative 
insight, we also present the phase diagram in a zero-hopping limit. In section HVl we present 
the phase diagram and the density profiles of 87 Rb- 41 K mixtures in a combined potential of 
a parabolic trap and an optical lattice. We will show that the phase separated 87 Rb- 41 K 
mixtures distribute a nesting structure. 



II. MEAN-FIELD APPROXIMATION FOR THE BOSE-HUBBARD HAMILTO- 
NIAN 

In this section, we describe the mean-field approximation for the Bose-Bose mixture in 
the absence of a parabolic potential. We consider two-component bosonic-atom mixtures 
trapped in an optical lattice at zero temperature. We use the following Bose-Hubbard 
Hamiltonian for the Bose-Bose mixtures; 



-t a y^filjbgj + b ] aj b ai ) + ^2 ^i(n ai - 1) - /i a ^n c 



+ ^12 J^li^i, (1) 

i 

where b^ ai (b a i) is the creation (annihilation) operator of the a{= 1 or 2) component at site 
i, h ai = b^bai is the number operator, and denotes the sum over the nearest neighbor 
sites. t a , U a and fi a are the hopping amplitude, the on-site intra-species interaction and the 
chemical potential respectively, and U12 is the on-site inter-species interaction. Throughout 
this paper, we consider the repulsive on-site inter-species interaction, i.e. TJ\i > 0. 
Using a mean- field approximation we decouple the hopping term as 

biA^(tyb a3 +bUk,)-(bl)(b aj ) 

= 4> a b a j + b ai (j) a - 4> 2 a , (2) 

where <p a = (S^) = (b a i) is the superfluid order parameter of the a component. Within this 
approximation, the mean-field Hamiltonian can be written as a sum over single-site terms, 
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H ~ Hmf = Y^ii H-h where the single-site Hamiltonian is given by 



-zt a (f) a (P ai + b ai ) + ^-h ai {h ai - 1) - fi a n ai + zt a (p 2 a 



+ U 12 h u h 2i , (3) 

where z is the coordination number. 

We can calculate the densities and the superfluid order parameters using the mean-field 
Hamiltonian ([3]). We employ the same basic method used by Lu and Xaing [4|, adapted 
to the case of Bose-Bose mixtures. The ground energy Eq and ground state wavefunction 
|^o) of Bose-Bose mixtures can be obtained by diagonalizing the Hamiltonian Hi in the 
occupation number basis {|ni, 7^2)} truncated at finite values na and n t2 . In the present 
work, we set n t \ = n t2 = n t and use a sufficiently large n t so that the results do not depend 
on n t . For a given U a , Ui 2 , t a , and /i a , the superfluid order parameters can be determined 
by minimizing Eq. The region with nonzero <p a is identified as the superfluid phase while 
the region with <p a = is identified as the Mott-insulator phase. After determining a , we 
obtained the superfluid densities p sa and average densities p a as 

Psa = 4%, (4) 
p a = (n a >. (5) 



III. PHASE-SEPARATION IN MI PHASE 



Before presenting detailed results of the mean-field theory, we clarify the condition for 
the phase separation. Here, we consider a zero-hopping limit t x (t 2 ) — > in order to obtain 
a qualitative understanding of the phase separation. We denote the energy of 1 711,712) per 
site by E ni ^ n2 , and we refer to this MI state as the (711,712) state. An explicit expression for 
this energy is given by 

E n u n 2 = y7ii(m - 1) + yni(n 2 - 1) + U 12 n 1 n 2 - Hxrii - fi 2 n 2 . (6) 

For the (ni,n 2 ) state to be the ground state, one must satisfy E ni ^ n2 < E ni ± ljTl2 , E niim ±i. 
This leads to the following relations between the chemical potentials and the number of a 
components; 

Ui(ni - 1) + Ui 2 n 2 < fii < U\n\ + U 12 n 2 , 

(7) 

U 2 (n 2 - 1) + U\ 2 n x < p 2 < U 2 n 2 + U\ 2 n x . 
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If /j, a < 0, there are no particles of the a component. 

We now consider the dependence on Ui 2 . In the limit U\% — > 0, the occupation numbers 
of two components should be independently determined by Eq. (|7|). On the other hand, in 
the opposite limit U± 2 —> oo, two components cannot coexist at the same site, so that the 
mixtures are phase-separated, i.e. the ground state is either the (ni,0) state or the (0,n 2 ) 
state. In this limit, the phase boundary between two phase-separated states is analytically 
given by: 



ni U 2 U\n x 

/i 2 = —Vi + -7r{n 2 - 1 - (m - I, 

n 2 2 2 n 2 



(8) 



where n\ and n 2 must be satisfied Eq. ([7]). 

We now derive the condition for the coexisting state. With decreasing magnitude of 
U\ 2 from U\ 2 — > oo limit, the coexisting state (jii, 1) or (l,n 2 ) first appears. In order for 
the (ni, 1) state to appear between the (rii,0) state and the (0,71.2) state, one must satisfy 
E nit i < E nit0 . Therefore, from Eq. ((6]) we obtain 



fj, 2 > U 12 ni. 

The intersection of fi 2 = U\ 2 ni and Eq.([H]) is given by 

/il = U 12 n 2 + -^-(ni - 1) - ^- — (n 2 - 1). 
2 2 ni 

Using Eq. (J7J), we find the condition for the appearance of the (rii, 1) state to be: 

1 



(9) 



(10) 



t/12 < 



n x n 2 



— ni(ni + 1) + — n 2 {n 2 - \, 



(11) 



From a similar calculation, the condition for the appearance of the (1, 112) state between the 
(ni, 0) state and the (0, n 2 ) state is given by 

1 



U 12 < 



— ni(ni - 1) + — n 2 {n 2 + l i 



nin 2 

These results are consistent with those of previous studies 
as in 



(12) 
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16|. If we set U x = U 2 = U, 



16|, the condition for the coexisting state is given as U\ 2 < U. We note that the 



condition for the coexisting state in the MI phase depends on the on-site interactions as well 
as the occupation numbers. 
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IV. NUMERICAL RESULTS 



In this section, we present numerical results of the mean-field calculation, considering 
the 87 Rb- 41 K mixture used in the experiment of Ref. 10] . We assume that a = 1 denotes 
87 Rb and a = 2 denotes 41 K. We calculate the densities and the superfluid order parameters 
using the method explained in Sec. [Til We estimate the on-site interactions and hopping 
amplitudes using the experimental data given in [loj]: lattice laser wavelength = 1064nm, 
scattering length agb = 99ao, clk = 65ao, am-K = 169ao, where ao is the Bohr radius. These 
results indicate that the on-site interactions are U\ : U2 '■ Un — 1 : 0-34 : 1.85, the hopping 
amplitudes are t\ : t 2 = 1 : 9.39, and U\ = AOzti. 



A. No hopping limit 

As an illustration, we consider the ground state in the no hopping limit, setting t a = 0. 
In Fig. [TJ we plot the phase boundaries of the ^1/^/1-/^2/^2 plane separating the regions 
with different occupation numbers. Figure DJa) shows the phase diagram obtained for the 
parameters described above. Since the on-site interactions do not satisfy the condition 
for coexistence, the mixture is always phase separated. Note that the slope of the phase 
boundary is different for different occupation numbers. The variation of the slope of the 
phase boundary is due to the difference of the on-site intra-species interactions U\ and {7 2 - 

We next consider the effects of changing the value of U\2- Experimentally, the magnitude 
of U12 can be controlled by the Feshbach resonance In Fig. Hfb), we plot the phase 

boundary for a weaker inter-species interaction U±2 = 0.6U±. In this case, the condition 
for coexistence (lllj) is satisfied, however, the condition (Tl2l) is not satisfied. Contrary to 
Fig. [T](a), the coexisting phase appears in the region enclosed by the solid line in Fig. QJb). 
However, most regions are still in the separated phases (ni,0) state or the (0,n 2 ) state. 



B. Finite hopping case (t a 7^ 0) 

We now consider the phase diagram for finite hopping amplitudes. In particular, we set 
the parameters for 87 Rb to be the MI phase and 41 K to be the SF phase. In Fig. [2} we plot 
the phase diagram of the 87 Rb- 41 K mixture at finite hopping. We find that the ground-state 
remains phase-separated, and the 41 K(SF) region becomes larger while the (ni,0) region 
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FIG. 1: The phase boundary of (ni,rt2) states in the zero-hopping limit with (a) U± : U 2 '■ U\ 2 = 
1 : 0.34 : 1.85, {h)Ui : U 2 : U 12 = 1 : 0.34 : 0.6. The solid line shows the boundary between 87 Rb 
and 41 K. The dotted lines show the boundary between the same species with different occupation 
numbers. 

becomes smaller than the t a — case. Since the SF phase appears in the ground-state, 
the phase boundary between the (ni, 0) state and the 41 K(SF) state now becomes a smooth 
curve in contrast to Fig. Ufa). 




FIG. 2: Phase diagram of Rb- K mixtures in the ^i/Ui-fi 2 /U 2 plane. The solid line represents 
the phase boundary of the {n\, 0) state and the (0, n 2 ) state. The dashed lines represent the phase 
boundary of the SF phase and the MI phase. 
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C. Effect of parabolic potential 



We now consider the effect of a parabolic potential. We will show that an interesting 
feature arises in the density profiles of 87 Rb- 41 K mixtures in a combined potential. We fix 
the number of 87 Rb and control the number of 41 K by changing the chemical potential \x-l- In 
order to include the effect of the trap potential, we apply a local density approximation (see 
Ref. Q). That is, we set the local chemical potential for each species as /^(ri) = fj? a — e Q (rj), 
where [x° a is the chemical potential at the center of a parabolic potential and e Q (rj) = ^^r? 
is the parabolic trap contribution for each species. It should be noted that we consider the 
difference in the masses of the two species. 

The variation in the density profiles is plotted in Figs. |3] (a)-(d). In Fig. |3]^a), we plot the 
density of a pure 87 Rb gas, i.e. 41 K atoms are absent. The density of 87 Rb in a parabolic 
potential exhibits the so-called wedding-cake profile, with a plateau corresponding to the 
MI (jii = 1) domain at the center of a parabolic potential and the SF state at the outside 
of the MI domain. Figure Efb) shows the density profiles of the 87 Rb and 41 K mixtures 
for fi^/zti = 16. Adding the 41 K component under the conditions for the phase separation 
f lTTT) and ( TT2l) . we see that 87 Rb exists in the center of a parabolic trap, which is surrounded 
by the 41 K component. An analogous structure is also found in the density profile of the 



BEC mixture without the optical lattice in a parabolic trap [18[. We note that the 41 K 
component in the outer region pushes the 87 Rb component. Thus, the Mott domain of the 
87 Rb component at the center of the trap, seen in Fig. Fig:DP (a), is turned into the SF 
region. This result is consistent with the results obtained using the Gutzwiller approach 
131 ] . We find that the Mott core changes to the superfluid at the center of a parabolic trap 
and 41 K exist outside the 87 Rb component. 

Upon further increasing the population of the 41 K component, the 87 Rb component of 
the SF state at the center of the trap is compressed, and thus the ri\ = 2 87 Rb Mott plateau 
appears at the center of the parabolic potential. However, not all the 87 Rb components 
change to the ri\ — 2 MI state and the n\ = 1 MI state remains. Moreover, the 41 K 
components penetrate between the n\ — 2 and n\ = 1 87 Rb Mott plateaus. Therefore, as 
can be seen in Fig. [3](c), the density profile becomes a nesting structure with alternately 
arranged 87 Rb, and 41 K from the center to the outside of a combined trap. This nesting 
structure is a peculiar feature of a Bose-Bose mixture trapped in an optical lattice. Upon 
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further increasing the population of the 41 K component, the ri\ — 1 87 Rb Mott plateau 
decreases and the ri\ = 1 87 Rb Mott plateau increases. In Fig. [3](d) , we plot the results for 
\i\/zt\ = 28. We find that the nesting structure seen in Fig. [3](c) disappears. As in Fig. E](b), 
all of the 41 K exists outside the 87 Rb component again. 
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FIG. 3: Density distribution of 87 Rb- 41 K mixtures in a combined potential. {a)ii^/zt\ = 0, 
(b)^2/zti = 16, (c)/!®/ zt\ = 25, and (cVjfi^/zti = 28. The solid line is the occupation number 
of 87 Rb, and the dashed line is the occupation number of 41 K. 



We can understand the nature of these density structures from Fig. HI In the local 
density approximation, the spatial dependence comes from the local chemical potential. 
The starting point of the trajectories of the chemical potentials for 87 Rb varies with 
changes in the chemical potential for 41 K. In the case of Figs. EJ^b) and (d), the trajectory of 
the chemical potential passes through the phase boundary of the (n\, 0) state and the (0, 722) 
state, as represented by the dotted and dash-dotted lines in Fig. HI On the other hand, in 
the case of the nesting structure, the chemical potentials start from the (n l7 0) region, pass 
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through the (0, 712) region, and enter the (774,0) region again, as shown by the dashed line 
in Fig. HI By the addition of more 41 K atoms, the trajectory changes to the dash-dotted 
line. The difference in the structure of the density distribution is caused by the difference 
in trajectories. 




FIG. 4: Phase diagram of 87 Rb- 41 K mixtures in the ^1/^1-^2/^2 plane. The solid line is the 
phase boundary of the (nx,0) state and the (0,77-2) state. We plot the trajectories of the chemical 
potential for fx^/zti = 16 (dotted line), fi^/zti = 25 (dashed line), and fj^/zti = 28 (dash-dotted 
line) respectively. 



V. CONCLUSIONS 

In this study, we investigated the properties of two-species Bose-Bose mixtures in an 
optical lattice using a mean-field approximation. We derived the analogous expression for 
the phase boundary between the (ni, 0) state and the (0, n 2 ) state by comparing the ground 
energy in the t a = and U\2 — > 00 limit. Using this phase boundary, we clarified the 
conditions for the spatial phase separation of Bose-Bose mixtures when both components are 
in the MI phase. These results depend on the on-site interactions and occupation numbers. 

We also studied the density distribution of 87 Rb- 41 K mixtures in a combined potential 
of a parabolic trap and an optical lattice using local density approximation. We estimated 



on-site interactions and hopping amplitudes using experimental data [IOJ, and calculated the 
densities and superfluid order parameters. We fixed the number of 87 Rb atoms and changed 
the number of 41 K atoms. We found that phase separation occurs even for 87 Rb- 41 K mixtures 
with 87 Rb as the SF phase and 41 K as the MI phase. We clarified that 87 Rb is localized at the 
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center of a parabolic trap and 41 K encircles the 8 Rb component. For a special case, we found 
a nesting structure, as shown in Fig. |3jc). These results are understood by considering the 
phase diagram in the ^1/^/1-/^2/^2 plane of the Bose-Bose mixtures given in Fig. [2J These 
results should be observed experimentally. 
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